Characteristic temperatures and spectral appearance of 
ULX disks 
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Abstract. A standard disk around an accreting black 
hole may become effectively optically-thin and scattering 
dominated in the inner region, for high accretion rates (as 
already predicted by the Shakura-Sunyaev model) . Radia- 
tive emission from that region is less efficient than black- 
body emission, leading to an increase of the colour temper- 
ature in the inner region, by an order of magnitude above 
the effective temperature. We show that the integrated 
spectrum has a power-law-like shape in the ~ 1-5 keV 
band, with a soft excess at lower energies and a downward 
curvature or break at higher energies, in agreement with 
the observed spectra of many ULXs. This scenario offers 
a physical explanation for the phenomenological "dense 
corona", successfully used to fit those spectra. It also re- 
lates their fit parameters to physical properties of the ac- 
creting black hole, such as mass and accretion rate. 



1. Introduction 

Fitting colour temperatures to the X-ray spectra of ac- 
creting black holes (BHs) has traditionally been a useful, 
indirect method to estimate BH masses, typically within a 
factor of 2. This method relies on the assumption that the 
accretion disk extends down to the innermost stable cir- 
cular orbit and is optically-thick, emitting a multicolour- 
blackbody spectrum (Shakura & Sunyaev 1973; Frank, 
King & Rainc 2002). If so, the fitted peak temperature and 
luminosity give a characteristic inner-disk area and hence 
BH mass. These conditions are usually verified in Galactic 
BHs in their high-soft state (0.01 ^ ix/^Edd < 1). 

For ultraluminous X-ray sources (ULXs), the scenario 
that the disk is optically thick and reaches the inner- 
most stable circular orbit could be a plausible initial 
assumption, in the absence of dynamic mass determi- 
nations. From a characteristic feature in the soft spec- 
trum ("soft excess"), peak disk temperatures « 0.15- 
0.20 keV have been inferred and used to suggest masses 
~ 1000M Q (Miller, Fabian & Miller 2004). A strong ob- 
jection against this argument is that most of the X-ray 
flux in ULXs comes out in a harder, power-law-like com- 
ponent, more consistent with inverse-Compton scattering 



(Stobbart, Roberts & Wilms 2006; Winter, Mushotzky & 
Reynolds 2006). This suggests that the soft thermal com- 
ponent only traces emission from the outer, optically-thick 
disk, at radii R ^ 10i?isco; if so, its temperature and 
luminosity cannot be directly used to estimate the BH 
mass. An additional parameter is needed to model the 
location of the transition radius between the outer, stan- 
dard disk and the inner, scattering-dominated region. In 
general, based on a variety of physical arguments, we ex- 
pect that the transition radius R c oc TTr -Rrsco f° r fh> 1, 
where m = M/Mgad an d /3 ~ 1. 

2. Alternative physical interpretations 

A comparison of X-ray spectral and timing properties of 
ULXs and of Galactic BHs, in particular XTE J1550-564, 
has suggested a possible physical interpretation in terms 
of a standard disk plus Comptonizing corona (Goad et 
al. 2006; Done & Kubota 2006). Unlike Galactic BHs, 
however, many ULXs have a second characteristic spectral 
feature at energies ~ 5-10 keV. There, the spectrum is not 
well approximated by a power-law: it shows a downward 
curvature or break (Stobbart et al. 2006; Roberts 2007). 
In the disk-corona framework, this can be interpreted as 
emission from a corona optically-thick to electron scatter- 
ing (r os ~ 10) and with a much lower temperature than 
is seen in Galactic BHs (kT ~ a few keV, compared with 
kT ~ 100 keV in Galactic BHs). In summary, the high- 
energy break would be related to the coronal temperature, 
and the soft component would come from the standard 
disk sticking out outside the corona, at large radii. 

A possible difficulty of this model is that, in order 
to be optically thick, the corona must also be extremely 
dense. Assuming a characteristic radial size of ~ 5000 
km, vertical size ~ 1000 km, r cs ~ 10 corresponds to 
n e ~ 10 17 cm -3 . This is many orders of magnitude denser 
than typical coronae in Galactic X-ray binaries. More im- 
portantly, one still has to explain how gravitational power 
is transferred from the underlying disk to the optically- 
thick corona. For proton densities n p rs n e ~ 10 17 cm -3 , 
the thermal energy stored in such a corona is only ~ 10 35 
erg, to be compared with a power ~ 10 40 erg s _1 per- 
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sistently transferred to the outgoing radiation via inverse 
Compton scattering. If gravitational power is released in 
the disk, a very efficient and stable energy transfer mech- 
anism is required to keep the corona energized. 

The same phenomenological spectral features have also 
been explained in a radically different way. The high- 
energy curvature and/or break have been interpreted as a 
very hot inner disk, with colour temperatures kT ln w 1.5— 
3 keV: e.g., Stobbart et al. 2006). A possible reason for 
such hot disks may be that the BH is rapidly spinning. 
Or, the disk may be non-standard, (e.g., slim disk mod- 
els: Watarai, Mizuno & Mineshige 2001). In this sce- 
nario, ULXs are a natural extension of Galactic BHs, with 
higher accretion rates and hotter disks (Fig. 3 in Stob- 
bart et al. 2006). The downside of this model is that it 
does not provide a natural explanation for the soft X-ray 
component. One possibility is that the soft excess comes 
from a downscattering outflow (King & Pounds 2003) or 
from relativistically-blurred emission lines in the reflection 
spectrum of a photoionized disk (Crummy et al. 2006). 

3. Modified blackbody spectra 

In most ULXs, both scenarios (cool disk with a hotter 
Comptonised component, or a hot disk with a soft excess) 
are consistent with the observed spectra. Before we can 
use the peak disk temperature as a mass indicator, we 
need to determine whether the disk is associated with the 
soft or the hard components. Here, we suggest that both 
interpretations are in a sense correct, and follow simply 
from the standard disk model (Shakura & Sunyaev 1973). 
The thermal emission spectrum of an accretion disk 
is successfully approximated with a multicolor blackbody 
model, provided that two assumptions are verified: the 
disk is optically-thick in the vertical direction, and all the 
gravitational power dissipated during the infall is radiated. 
Then (see, e.g., Frank et al. 2002 for details), 
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where D(R) is the power dissipated per unit disk area, 
d is the distance to the source, and F° bs is the observed 
flux. Before substituting Eqs. 1 and 2 into Eq. 3, one has 
to make sure that the disk is optically thick at all R. The 
effective optical depth (Rybicki & Lightman 1979) 
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where H(R) is the disk thickness, takes into account both 
the (frequency-dependent) free-free absorption coefficient 

a5 « 3.7 x 10 8 T- 1 / 2 Z 2 „ e r» i i/-3(i - e- hv ' kT )g & cm" 1 (5) 



(where the Gaunt factor g 
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1), and the electron scat- 
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At large radii, for R £ 7.2 x 10 9 (m) 2 ^(M/M Q ) cm 
(e.g., Frank et al. 2002) a standard disk is gas-pressure 
dominated, optically thick (r° ff ^> 1), and the Rosseland 
mean absorption opacity Tp 3> r cs . In this zone, I v = 
S v = B v and the spectrum emitted by each annulus is a 
blackbody. 

Moving towards smaller radii, there is a region where 
the disk is still gas-pressure dominated and optically thick 
(r° ff 3> 1) but the opacity is dominated by electron scat- 
tering (t° s ^> t5 ^> 1). Deep inside the disk, the radiation 
is still in thermal equilibrium with the matter (S v = B u ), 
but the emerging spectrum from the disk surface is mod- 
ified by scattering, such that the outgoing intensity 
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(Rybicki & Lightman 1979). From Eq. (7) one can directly 
calculate the integrated disk spectrum (see Section 4). It 
is often useful to introduce an approximation of Eq. (7) 
that does not contain a frequency dependence, by using 
the Rossland mean of the absorption coefficient 
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where the frequency-averaged Gaunt factor j^rjI. Then, 
the total emitted flux from a unit area on the disk 
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that is, for a given temperature, the disk is slightly less 
efficient than a blackbody at radiating the dissipated grav- 
itational power. If we impose that the disk has to radiate 
the same power, the colour temperature at each radius 
must increase with respect to the effective temperature: 



T'(il)«0.81(a"7a!) 1 / 8 T eff (/e), 



(10) 



where 1 < (a cs /a^) 1,ls 5s 2 as we will show later. 

At smaller radii, R < 1.1 x 10 8 a 2/21 m 16 / 21 
(3) (M/M ) 13 / 21 cm (Frank et al. 2002), radiation pressure 
dominates over gas pressure. In this region, the disk may 
become optically thin to true absorption (r* <C 1), but 
is always optically thick to scattering (t° s 3> 1) and most 
importantly, is still effectively optically thick (r° ff 3> 1). 
As long as r^ ff 3> 1, Eqs. 7 and 9 are still valid, and the 
emerging spectrum is still a modified blackbody. Explic- 
itly, for a pure Hydrogen gas, from Eqs. 5 and 6: 



1.8 x lO-^T^VVp - e- hv ' kT )-\ (11) 

(H/ cos i) J a cs a^ w 1.6 x 10~ 8 (H/cosi) 
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Applying standard-disk solutions for the radiation- 
pressure dominated region (Shakura & Sunyaev 1973; 
Frank et al. 2002), we can approximate 

T{R) w 2.3 x lO 7 ^- 1 / 4 ™- 1 /^- 3 / 4 K, (13) 

H(R) w 2.5xl0 6 ram[l-r" 1/2 | cm, (14) 

n e (R) S3 6.1 x lO^o^m^mrV 3 / 2 
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where we have defined m = M/Mq, r = R/Risco- Insert- 
ing Eqs. 13-15 into Eqs. 11 and 12, we obtain: 
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For typical Galactic BHs in a high state, the innermost 
part of the disk falls into the radiation-pressure, electron- 
scattering-dominated, effectively-thick regime. For M — 
10M Q , R = 2i?isco, M = -^Eddi viscous parameter 
a = 0.1, hv = 1 keV (so that 1 - exp (-hv/kT) ~ 1 
across the 0.3-10 keV band), we see from Eq. 16 that 
the colour correction ,f u {R) = .81(a cs / 'af,) 1 / 8 S3 2, only 
weakly dependent on frequency (within the X-ray band), 
accretion rate, BH mass and radius. In fact, 1.5 ^ /„ Ss 2.5 
throughout the X-ray emitting regions, for the character- 
istic range of physical parameters typical of the high/soft 
state in Galactic BHs. That is why f u (R) can be approx- 
imated with a constant "hardening factor" / S3 1.7-2.5 
(Shimura & Takahara 1995; Shafee et al. 2006; Shradcr & 
Titarchuk 2003). 

The previous results depend critically on the condi- 
tion T° ff <; 1. From Eq. 17, using the same set of param- 
eters suitable for Galactic BHs in a high state, we see 
that t° h S3 5 in the innermost region, at the Eddington 
accretion rate; at larger radii, t„ ^> 1. So, the optically- 
thick condition is always verified for Galactic BHs, but 
the X-ray emitting region gets close to becoming effec- 
tively thin. A small increase in viscosity, BH mass and/or 
accretion rate above Eddington will lead to the formation 
of a radiation-pressure dominated, effectively-thin region 
in the inner disk (r® Ss 1, with r^ <C 1 and r° s <; 1). 
If that happens, Eq. 7 is no longer applicable. Note that 
an increase of the accretion rate leads to a reduction of 
the effective optical depth in the radiation-dominated zone 
(perhaps counter- intuitively) . This is essentially because 
n e ex mT 2 . Conversely, in the gas-pressure dominated disk, 
an increase in m leads to an increase in optical depth. 

In the region where the disk becomes effectively thin, 
the specific intensity emitted by an annulus is 
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If all the flux is radiated from each annulus, the temper- 
ature must increase even more dramatically than in the 
optically-thick regime. The hardening factor is now 

f v {R) S3 (a^/cosz)" 174 

S3 21.0i/?£ v (1 - e- fc "/W)-i/4 (i/cos*) 1 / 4 
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Let us consider a hypotetical ULX with a mass m = 
100 and an accretion rate m <; 1. Then, the inner disk 
is effectively thin for "reasonable" choices of a, and the 
transition radius between thin and thick regions is 

r c £ 12.1a 23 / 39 m 7 / 39 m 32 / 39 (l - r- 1 ' 2 )^'™ (cos*) 16 / 39 
8 . (20) 



a few 



For r <; r c , the spectrum is a disk-blackbody, slightly mod- 
ified by a hardening factor S3 constant over the frequencies 
and radii relevant to the X-ray spectral observations, and 
most likely Ss 2. For r Ss r c , the hardening factor increases 
sharply towards the innermost stable circular orbit, reach- 
ing /ikeV ~ 4.3 to 3 ' 4 <~ 10 in the brightest region of the 
disk. Although the effective optical depth is < 1 in this 
region, the optical depth to scattering alone 
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4. Observational predictions 

If a disk develops an effectively-thin inner zone, how would 
its spectrum differ from a multicolour blackbody? To 
have a qualitative idea, let us consider for simplicity a 
frequency-averaged hardening factor f(R) > 1, such that 
I V {R) = f- i {R)B v [T{R)\ and T{R) = f(R)T cS (R). Then, 
the observed flux can be obtained as usual by integrating 
the emission over all disk annuli (Eq. 3): 
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Using for example the hardening factor plotted in Fig- 
ure 1 (top panel), where we have assumed that / = 1 for 
r > r c , and applying Eq. 22, we obtain the X-ray spec- 
trum plotted in the bottom panel of Figure 1. A detailed 
study of the emitted spectra is left to a paper currently in 
preparation. Here we simply want to point out the main 
qualitative result: the spectrum has a power-law-likc shape 
for energies 1.5 keV Ss E ^ 6 keV, a downward curva- 
ture and break above 6 keV, and a "soft excess" with a 
characteristic temperature T c ~ 0.2 keV. The latter com- 
ponent corresponds to the disk-blackbody emission from 
the outer disk (R > R c ), where / ~ 1. The high-energy 
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Energy (keV) 

Fig. 1. A hardening factor that increases sharply at small 
radii (top panel), as we suggest is the case in ULX disks, 
produces an emitted X-ray spectrum (bottom panel) with 
three characteristic features: a power-law-like component, 
a soft excess at low energies, and a curvature or break at 
high energies. 



above the observed X-ray band, respectively, creating an 
apparently "pure power-law" spectrum. 



5. Conclusions: standard disk or corona? 

We have shown that the accretion disk in BHs with 
M <; lOAf , M ~ a few MEdd is expected to develop 
a scattering-dominated, effectively-thin inner region (R ^ 
10 i?isco)> with a corresponding increse in the colour tem- 
perature, to compensate for the decreased emission effi- 
ciency. The predicted spectrum is very similar to what is 
observed from many ULXs at Lx ~ 10 40 erg s _1 . Phe- 
nomenologically, those spectra have been successfully fit- 
ted with an optically-thick disk-blackbody emission from 
the outer disk, plus a dense scattering corona in the in- 
ner region, "covering" the inner disk. We have shown here 
that the characteristic scattering opacity, size and tem- 
perature observationally inferred for those "coronae" are 
the same as predicted for a radiation-pressure dominated 
standard disk, for accretion rates and BH masses consis- 
tent with ULX parameters. So, we argue that there is no 
need to postulate an ad-hoc corona covering or replacing 
the inner disk, ft is the standard disk itself that becomes 
a hot, scattering-dominated environment at small radii, 
and produces the observed spectra. The advantage of see- 
ing it as the inner disk rather than a corona is that we do 
not have to worry about how energy is transferred from 
one to the other, or how their interface looks like. More- 
over, in the scattering disk scenario, the peak tempera- 
ture, optical depth and emitted flux are not just empirical 
parameters: they are predictable functions of M, M, a, 
known from standard-disk models. This provides the tools 
to determine the physical properties of the accreting BHs 
in ULXs from their fitted spectral parameters, and relate 
their spectral evolution to changes in the accretion rate. 
In this scenario, both the cool (T ~ 0.2 keV) and the hot 
component (T <; 5 keV) come from the disk. 



break corresponds to the maximum colour temperature 
in the effectively-thin, scattering-dominated inner disk: 
T rs fT e g m 6 T c ff- This result is already well known 
(Shimura & Takahara 1995; Fig. 12 in Shakura & Sunyaev 
1973), but its implications for ULX studies are not well 
explored. In between the two thermal features, the slope 
of the power-law-like section of the spectrum depends on 
how steeply the hardening factor increases towards the 
centre. A slightly steeper slope is obtained when the fre- 
quency dependence in the hardening function /„ is also 
taken into account. Additional steepening of the spec- 
trum results from non-radiative contributions (outflows, 
Poynting flux) to the cooling flux from the inner region. 
Increasing the range of values of the hardening factor (for 
example, by increasing the viscosity or the accretion rate) 
moves the soft excess and high-energy break below and 
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